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We study charmonium correlators in pseudoscalar and vector channels at finite temperature using 
lattice QCD simulation in the quenched approximation. Anisotropic lattices are used in order to 
have sufficient numbers of degrees of freedom in the Euclidean temporal direction. We focus on the 
low energy structure of the spectral function, corresponding to the ground state in the hadron phase, 
by applying the smearing technique to enhance the contribution to the correlator from this region. 

We employ two analysis procedures: the maximum entropy method (MEM) for the extraction of the 
spectral function without assuming a specific form, to estimate the shape of the spectral function, 
and the standard fit analysis using typical forms in accordance with the result of MEM, for a 
more quantitative evaluation. To verify the applicability of the procedures, we first analyze the 
smeared correlators as well as the point correlators at zero temperature. We find that by shortening 
the t-interval used for the analysis (a situation inevitable at T > 0) the reliability of MEM for 
point correlators is lost, while it subsists for smeared correlators. Then the smeared correlators at 
T ~ O.QTc and l.lTc are analyzed. At T ~ O.QTc, the spectral function exhibits a strong peak, well 
approximated by a delta function corresponding to the ground state with almost the same mass as 
at T = 0. At T ~ l.lTc, we find that the strong peak structure still persists at almost the same 
place as below Tc, but with a finite width of a few hundred MeV. This result indicates that the 
correlators possess a nontrivial structure even in the deconfined phase. 

PACS numbers: 12.38.Gc, 12.38.Mh 


I. INTRODUCTION 

The charmonium systems have been payed much at¬ 
tention as a signal of the QCD phase transition. Due to 
a change of the interquark potential by thermal effects, 
the J/ip mass is expected to decrease when approaching 
the phase transition P|. Above Tc, the screening of the 
interquark potential may dissolve the charmonium states, 
and the resulting J/ifi suppression has been regarded as 
one of the most important ^nals for detecting the for¬ 
mation of the plasma state [illli 

In principle, lattice QCD simulations can provide in¬ 
formation on such excitation modes based on a nonper- 
turbative and model independent framework, since the 
correlation functions in the Euclidean temporal direction 
measured on a lattice are related to the real time re¬ 
tarded and advanced Green functions by analytic contin¬ 
uation from a single spectral function |^ |^ . In practice, 
however, the extraction of reliable information from nu¬ 
merical data for correlators becomes increasingly difficult 
as the temperature increases. There are two restrictions 
on the correlators at finite temperature: the numbers of 
available degrees of freedom, and the physical extent in 
the temporal direction. The former can be improved by 
the use of anisotropic lattices, on which the teinporal lat¬ 
tice spacings are finer than the spatial ones [a H 111 ■ The 
latter restriction is, however, physically inevitable and 
may spoil standard procedures applicable at zero tem¬ 
perature. 

Recently, there has been technical progress in the di¬ 
rect extraction of the spectral function from lattice data 


of correlators UM The key role is played by the maxi¬ 
mum entropy method (MEM), combined with the singu¬ 
lar value decomposition for constructing the functional 
bases in the space of spectral functions. At zero temper¬ 
ature this technique has reproduced ^ectra consistent 
with the standard fit analyses . Appli¬ 

cation to systems at finite temperature is, however, not 
straightforward because of the two restrictions mentioned 
above [Ifil. In particular, shortness of temporal extent, 
1/T, makes it difficult to extract the correct low energy 
structure of the spectral function. Therefore, one at least 
needs to verify at T = 0 that the method produces results 
which are stable against variation of the Tinterval. There 
have already been a number of app lications of MEM for 
finite temperature Elllllll Il9l |. but as long as such 
tests are not systematically performed these results may 
contain uncontrolled uncertainties. In fact, we will show 
in this paper that MEM applied to point-point correla¬ 
tors at T = 0 with shortened t-interval fails to produce 
the correct spectral functions. We therefore apply the 
smearing technique to enhance the contribution of the 
low energy part to the correlators Him For the 
smeared correlators at T = 0, MEM works, at least qual¬ 
itatively, also with a shorter Tinterval. 

In this paper, charmonium correlators are investigated 
using lattice QCD simulations in the quenched approx¬ 
imation. We focus on properties of the ground state, 
more generally the low energy structure of the correla¬ 
tor, rather than the whole spectral function. This paper 
mainly deals with the following two subjects: 

(1) A technical investigation of procedures for extrac¬ 
tion of reliable information from the numerical data of 
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correlators in the temporal direction. In addition to 
MEM, we apply the standard fit anal ysis assuming 
a peak structure for the spectral function [^. The lat¬ 
ter approach complements the former, since MEM with 
small numbers of data points appears not work beyond 
the level of qualitative estimate for the shape of spectral 
function. 

(2) A study of temperature effects on the correlators 
near to the deconfining phase transition. Below Tc, the 
question is whether the charmonium mass shifts, as ex¬ 
pected from the potential model approach Q. Above Tc, 
previous studies have indicated that nontrivial structures 
may survive HIS The possibility of the existence of col¬ 
lective modes in the plasma phase is therefore carefully 
examined. 

We use quenched anisotropic lattices with the renor¬ 
malized anisotropy ^ jar = 4 and the spatial lattice 
cutoff ~ 2 GeV. The phase transition occurs near 
Nt = 28, and hence sufficiently many points are available 
in the vicinity of Tc. At T cs 0 we try to determine the 
conditions for a reliable extraction of the spectral func¬ 
tion. We find that for the point correlators, reliability 
of MEM with less than 24 t-points (meaning a physical 
Grange of less than 0(0.5fm) ~ 0.7Tc) is not guaranteed 
with the present level of statistics. Therefore only the 
smeared correlators are analyzed at finite temperature, 
T ~ 0.9Tc and 1.1 T^. 

Here we comment on the smearing technique. Al¬ 
though the correlator of spatially extended operators 
has no direct relation to the physically observable dilep¬ 
ton production rate, the properties of collective excita¬ 
tions such as mass and width are unchanged and can be 
probed more efficiently. To the extent that these excita¬ 
tions drive the dilepton rate we can at least say some¬ 
thing about their position and shape, if not about their 
strength. A disadvantage is the possibility of detecting a 
fake peak, which may be produced by the smearing even 
in the case of free quarks [T^. To distinguish such an 
artifact from a genuine physical peak we use besides the 
smearing function based on the wave function obtained 
at T = 0 also a narrower one of essentially half-width. 
We speak thereby of “smeared” and of “half-smeared” 
correlators (because the energy region enhanced by this 
narrower smearing function is wider, while the high fre¬ 
quency part of 0 ( 1 / 0 ,-) is still sufficiently suppressed). 
Comparing the results of these two types of smeared cor¬ 
relators, one can argue whether the observed peak struc¬ 
ture is an artifact of free quarks or physical indication of 
a collective mode. 

This paper is organized as follows. In Sec. im we re¬ 
call some fundamental properties of the spectral func¬ 
tion. For later uses in the analysis, presumable 

forms of the spectral function for a collective mode are 
discussed. Section Hill reviews anisotropic lattice actions 
and describes the parameters used in this paper. We 
also show the spectral function of the correlator of a free 
quark-antiquark pair, for comparison with the correlators 
from the Monte Carlo simulation. In Sec. CYI our anal¬ 


ysis procedures are explained. Section El describes the 
setup of Monte Carlo simulation. The following three 
sections show the result of the numerical simulation. In 
Section ED we analyze the correlators at zero tempera¬ 
ture, as a preparation for the more involved situations 
at finite temperatures. Sections EID and Emi show the 
results at T cs 0.9Tc and l.ITc, respectively. In each 
case, we first estimate the presumable form of the spec¬ 
tral function with MEM, and then apply the fit anal¬ 
ysis for a more quantitative evaluation of its structure. 
Section is devoted to our conclusions on the technical 
and physical implications of the results. Our preliminary 
results have been reported in Ref. . 

II. SPECTRAL FUNCTION 

In this section (3 represents the inverse temperature 
and should not be confused with the coupling of gauge 
field which appears in later sections. For a finite lat¬ 
tice of temporal extent Nt we have (3 = 1/T = NtUr- 
The Euclidean time is represented by r in this section 
to make clear the distinction between the Euclidean and 
Minkowski formulations. In later sections we use also t 
as the Euclidean time, because no confusion is expected 
for the lattice QCD formulated in the Euclidean space. 

We consider correlators of the form: 

Cir)=Y,{O{x,r)OH0,0)}, ( 1 ) 

X 

where the operator 0 (x, r) is a quark bilinear, 

0{x,t) ='^(j){y)q{x + y,T)Tq{x,T). (2) 

V 

4x4 matrix T specifies the quantum number, such as, 
for example, T = 75 for pseudoscalar and 7 ^ (i = 1,2,3) 
for vector channels. is a smearing function defining 
the spatial extension of the bilinear operator. Since the 
smearing in general violates the gauge invariance of the 
correlator, one needs to fix the gauge or employ a gauge 
invariant form of the smearing function. We fix the con¬ 
figurations to the Coulomb gauge, since it is widely used 
and allows easy implementation of the smearing func¬ 
tion. Source and sink operators are always identical in 
this work. 

Point correlators correspond to a ultra local “smear¬ 
ing function” ^(y) = ia Eq- For the smeared 
correlators we use smearing function of the form 

cj){y) = exp(-a|y|^’), (3) 

where the parameters a and p are determined by fitting 
the wave function measured in the numerical simulation 
to the above form. The smearing technique described 
here was already applied to the problems of the char¬ 
monium correlators at T > 0 |^, as well as to the light 
meson correlators 0 - 
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For a while let us consider the continuum Euclidean 
space. The Matsubara Green function is represented in 
terms of the spectral function as 


G{Vn,P} 


r°° du p{uj,p) 

27r iVn-bj + it' 


( 4 ) 


The spectral function p(ijJ,p) is represented as 

' n,m 

x{2Tr)^S{uj -Em+ E„)S{p- + P„), (5) 


which is an odd real function due to the bosonic nature 
of the present correlator and positive for w > 0 . 

The retarded and advanced Green functions for real 
time are represented as 


G«(z.,p) = 

duj 


(6) 

i -00 27r 

V — UJ + ie 

g^(gp) = 

r°° duj 

p{uj,p) 

(7) 

i-oo 27r 

u — UJ — ie 


with the same spectral function, p{uj,p) 0]. The spectral 
function is then represented as 

p{i',p) = —2 ImG^(z^,p) = 2 lia.G^{v,'p). ( 8 ) 


Strictly speaking, the above integrations do not always 
converge, and one needs appropriate subtraction terms. 
On the other hand, the high frequency part of the spec¬ 
tral function is not practically significant in the numerical 
analysis, because of the strong suppression by K{t,uj), 
and of the existence of the lattice cutoff. We therefore 
do not consider these subtractions in this paper. We 
note that the smearing of the operator in Eq. © is per¬ 
formed only in the spatial directions, and therefore the 
analytic continuation of the Matsubara Green function 
to the real time Green functions is independent of this 
operation (but, of course, the spectral function depends 
on the operators, hence on smearing). 

In the following we are only concerned with the correla¬ 
tors projected on the zero momentum states and discard 
the argument p. Since the spectral function p(uj) is an 
odd function in uj, the correlator o is represented as 

pOC 7 

C('r)=y^ ^K{t,lv)p{u;). (9) 

where 

p-UJT -0U1+UJT 

■ ( 10 ) 

In the numerical simulation of lattice QCD one measures 
the left hand side of Eq. ©■ The question is how to solve 
the inverse problem to obtain p{uj) with limited data for 
C{t). The procedures to extract information about p{uj) 
are described in Sec. nvi 


Now we shall present several ansdtze for the spectral 
function which will be used for later analysis of the lattice 
data, and discuss their physical implications. A simple 
representation of a collective mode is the following form 
of the retarded Green function: 


G^{u) 


V - p + 


( 11 ) 


where p and ^{v) express the dispersion and the width 
of the mode. The corresponding spectral function is the 
well-known Breit-Wigner form. We shall neglect the fre¬ 
quency dependence of the width 7 for simplicity. Taking 
the oddness of the bosonic spectral function into account, 
p{lo) reads 


p(^) =___. (12) 

’ {lo- mf (w-km) 2 -H 72/4 ^ ’ 


If the quantum numbers of the operators used to repre¬ 
sent the collective mode are the correct ones, the physical 
properties associated with the mode, characterized by p 
and 7 , are independent of the particular operator. There¬ 
fore the smeared operators can be used to observe these 
quantities. However, the smearing does change the over¬ 
lap of the operator with the state, i?(a;). Therefore the 
effects depending on and those coming only from 

the peak structure (represented with p and 7 ) should be 
distinguished. If the change of is sufficiently gentle 
over the region of the observed peak, the uj dependence 
in R{u}) is negligible. 

Instead of Eq. at zero temperature one often uses 
the relativistic Breit-Wigner type form: 


p{uj) 


sign(a;)w^m 7 ( 4 i?/m) 
(aj2 — m2)2 + TO272 


(13) 


This form shows a similar behavior to Eq. m around 
the peak position. However, the behavior far from the 
peak is different. In particular c© and ca give differ¬ 
ent contribution to the integral m near vanishing to: In 
the vicinity of cj = 0, Eq. 03 behaves linearly in w, while 
Eq. (ICT behaves linearly in Since K{t,uj) behaves as 
{2 — Idoj)/l3oj, Eq. {13 gives a t-independent contribution 
to the integral © p)a . We note that there is no indica¬ 
tion of linear behavior in uj in the small uj region of the 
spectral functions obtained from our numerical data with 
MEM, as shown in Sec. IVHI and |VH1I Therefore, in the 
fit analysis, we use the form of Eq. da as ansatz for 
the spectral function to which the lattice data are fitted. 
While the use of Eq. da is valid at zero temperature, 
it should only be taken as a representative form of peak 
structure at T > 0 , and hence the parameters m and 7 
do not have a strict sense but represent only a convenient 
parameterization of the peak. 

In the limit of 7 ^ 0, both Eqs. da and da tend 
towards a delta function form, 


p{uj) = 2 'kR [J(a; — m) — 5{uj + m)]. (14) 


This form corresponds to the standard exponential fit 
of correlators used in the hadron spectroscopy at zero 
temperature. 
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III. ANISOTROPIC LATTICE QCD 


Instead of {ko-, Kr), it is convenient to use the parameters 


A. Anisotropic lattice actions 

For analysis of the correlators in the temporal direction 
at T > 0 the fine temporal resolution is very important 
to extract meaningful information Anisotropic 

lattices have become a powerful tool to achieve suffi¬ 
ciently high temporal cutoff while keeping total compu¬ 
tational size modest. In the following, we summarize the 
anisotropic lattice actions employed in this work. The 
parameters we use are described in the next subsection. 

For the gau ge fi eld, we adopt the standard Wilson pla- 
quette action 


JF = UrKr/UaUcr, (19) 

- =-2(7f- b 3r — 4) = 2(mo7F + 4), (20) 

K Ua^a 

where mg is the bare quark mass in temporal lattice 
units. K plays the same role as on isotropic lattices, and 
controls the bare quark mass. 7 f is the bare fermionic 
anisotropy. 

The lattice quark field q{x) in the action Ijf 7|l is related 
to the dimensionful field ip{x) as 

^{x) = K)K:{K)q{x), (21) 


3 r 


= / 37 gEE 


X i—1 


1 - -ReTrC/i4(a;) 


where we assume m n -C JC{k) is so-called KLM 

normalization factor |2 (t|| . 

JC{k) = [2 mt-Kt-( 1-b mo)]^/^ 


X i<j^l 


1 - ]-ReTTUij{x) 

, (15) 

1 -b 2k{'Jf - 4) 

L 3 J 


1 -b 2k(^f -b 3r — 4) 


-,1/2 


( 22 ) 


where 1/^;^ is a product of link variables along a plaquette 
in the plane. Here the parameters (3 and yc are bare 
coupling and bare anisotropy, respectively. 

For the quark field we use the 0(a) improved Wilson- 
type action [^ 1^1^ : 

Sq = '^q{x)K{x,y)q{y), (16) 

x,y 


K{x,y) = 6^y 

-Kr [(1-74)04 (x)< 5,^+4 -b (l-b74)04(a; - 4)5^_4 
3 

[ir-li)Uiix)S^^i^y + ir+-f^)U,{x-i)6^_iy 


3 

K(jCf ^ ^ ^ 
i=l 
3 

K(yCB ^ ^ ^ijFij{x')Sx^y- ( 47 ) 

i>j=l 


where Ko- and Kr are the spatial and temporal hopping 
parameters, respectively, r is the spatial Wilson param¬ 
eter, and ce, cb are the clover coefficients for the 0(a) 
improvement. We set r = 1/^ and ce, Cb to the tadpole- 
improved tree-level values, and vary the two parameters 
Ka and Kr to change the quark mass and the fermionic 
anisotropy. The tadpole improvement |2 !t| is performed 
by rescaling the link variables as Ui{x) Ui{x)/ua and 
04 (x) ^ U 4 {x)/ur, respectively, with the corresponding 
mean-field values of the spatial and temporal link vari¬ 
ables, Ucr and Ur- Then ce and cb with the choice r = 1/^ 
read as 


(18) 


We consider the tadpole-improved tree-level matching of 
the fields in continuum and lattice theories, and hence 
Zq is set to unity. 


B. Parameters 

We here describe the parameters used in this paper. 
First, we discuss the calibration of gauge and quark fields. 
In an interacting theory, the renormalized anisotropy 
^ = Uajur generally differs from the bare anisotropy pa¬ 
rameters 7(3 and y^. Therefore one needs to tune the 
latter such that 

Cg(7g,7f) = Cf(7g,7f) = C (23) 

holds, where and are the observed anisotropies 
defined through the gauge and fermionic observables, re¬ 
spectively. In quenched simulations, one can calibrate 
the gauge and fermion fields separately, first the former, 
and then the latter. 

For the calibration of the gaug:e field, we refer to 
an elaborated work by Klassen which uses Wilson 
loops in the spatial-spatial and spatial-temporal planes 
and is performed at 1% accuracy level. Making use 
of his relation of yc to (3 and we choose (3 = 6.10 
and 7G = 3.2108 corresponding to the renormalized 
aniso trop y ^ = 4. These values were also adopted in 
Ref. 13, and correspond to the spatial lattice cutoff 
a~^ = 2.030(13) GeV set by the hadronic radius r n 1^ . 

For the quark field, we use the result of Ref. in 
which the calibration was performed on the same lat¬ 
tice as this work using the meson dispersion relation in a 
quark mass region including the charm quark mass. Ac¬ 
cordingly, we adopt k = 0.112 and y^ = 4.00, namely 
vriq = 0.121, as the values corresponding to the charm 
quark mass. (In Ref. it was argued that the value of 


Ce = I/mctM^, 


Cb = l/ul- 
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^ E 

i.fc=i.2 ■- 


Vj{-p)ruk{p) 


1 2 


2E{p) 


S (a 


2i?(p)X24) 


where Uj(p) and Vj(j>) are the quark and antiquark 
spinors with j-th spin, and (/>(p) is the Fourier transform 
of </>(?/). The summation in p is taken over all modes 
on the finite lattice of the same size as in Monte Carlo 
simulation, 20^. With the present lattice quark action, 
the energy of free quark, E(j)), satisfies the dispersion 
relation 


cosh if = 1 + 


mo + ^P 


l + mo + 


(25) 


FIG. 1: The spectral functions for the pseudoscalar correla¬ 
tors composed of free quarks at T = 0. 


7f tuned for the massless quark is also applicable to the 
charm quark mass. For historical reason, however^ we 
use the result of the mass dependent tuning in Ref. |^ . 
The difference is actually small and does not cause any 
significant effect.) 

As smearing function we use the result for the wave 
function in the Coulomb gauge at T = 0. For the vector 
channel the fit to Eq. @ yields a = 0.2275(9) and p = 
1.258(5). These values are used for both the pseudoscalar 
and vector correlators. 

In addition, we also use a narrower smearing func¬ 
tion with a = 0.45 and the same p. The correlators 
smeared with this function are called half-smeared corre¬ 
lators, since the smearing function has about twice the 
slope as the main one and hence a smaller smearing effect 
is expected: the low energy part of the correlator is only 
partially enhanced, while the high energy part of 0(l/ai-) 
is sufficiently reduced. This intermediate smearing func¬ 
tion plays an important role at T > Tc in examining 
whether the observed peak structure of the spectral func¬ 
tion for the smeared correlator is an artifact of smearing 
or a genuine physical effect fSec. IVIlUl . 


C. Free quark case 

Since in the perturbative plasma regime we expect free 
quarks and gluons it is important to compare the result 
of Monte Carlo simulation with the case of free quarks 
on the lattice. The spectral function of the correlator m 
in the case of free quarks is obtained from the expression 
0 . In the case of vanishing total momentum, for positive 
frequency, 

p 


where pi = (1/7f) sinp^ and pi = 2sin(pi/2). We set the 
bare quark anisotropy 7 f = ^ = 4 and the bare quark 
mass Too = 0.121 such that the free quark has the same 
value as in the Monte Carlo simulation. 

Figure Q shows the spectral functions of point, 
smeared, and half-smeared correlators composed of free 
quarks in PS channel. The integral of these functions 
are normalized to unity. Since on the finite lattice the 
spectral function is a sum of delta functions, we repre¬ 
sent each contribution at momentum p with a Gaussian 
of width 0.005. As clearly shown in the figure, the smear¬ 
ing of the operator strongly enhances the low energy part 
of the spectral function, as compared to the point corre¬ 
lator. In the (fully) smeared correlator, several low mo¬ 
mentum states over the range of uj about 0.1 dominate 
the correlator. From this correlator, therefore, an analy¬ 
sis with poor energy resolution in lo might produce a peak 
with full width of order of 0.1 instead of a sum of distinct 
peaks corresponding to individual momentum states. As 
we will see in Sec. ITvTtI the fit analysis with the single 
Breit-Wigner type ansatz applied to the smeared corre¬ 
lator composed of free quarks produces a much narrower 
peak with the width of about 0.027 centered at w ~ 0.275. 
This is of the same order of magnitude as the width ex¬ 
tracted from the correlators in Monte Carlo simulation 
at T ~ l.lTc. Therefore this occasion must be examined 
carefully. 

It is for this purpose that the half-smeared correla¬ 
tor is introduced. As shown in Fig. ^ a slightly higher 
and wider frequency region contributes dominantly to the 
half-smeared correlator than to the smeared one. If the 
correlator is approximated by a single peak, the peak 
will be centered at larger energy and have wider width as 
compared to those of the smeared correlator. This behav¬ 
ior is verified in Sec. CUby applying the fit analysis to 
the correlators composed of free quarks. Therefore, com¬ 
paring the results for spectral functions for the smeared 
and half-smeared correlators can help judging whether 
the peak observed from the simulation data is a genuine 
excitation mode or an artifact brought by the smearing. 
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IV. ANALYSIS PROCEDURE 
A. Strategy of analysis 

One of the goals of this work is to investigate tech¬ 
niques to extract the relevant information of the spectral 
function from the correlators with limited numbers of the 
degrees of freedom. We treat an inverse problem repre¬ 
sented as 

poo 

C{t) = / dujK{t,u>)A{uj), (26) 

Jo 

where C{t) is the given lattice result for the correlator, 
and the spectral function Aiuj) is what we need to ob¬ 
tain. Hereafter we shall denote the spectral function re¬ 
constructed from the lattice data A{uj) instead of p{lo) 
(x(27r)“^). Unless stated otherwise, the variables are 
in temporal lattice units in the following. The kernel 
is given as 


combination provide a more reliable way to analyze the 
structure of spectral functions than taken independently. 
We thus propose a two-step procedure: we first apply 
MEM to the correlators, for a rough estimate of the shape 
of spectral function. Once a presumable form is known, 
we then examine this form using fit, and estimate the 
values of parameters such as the mass and width of a 
peak structure. 

B. Maximum entropy method 

Our MEM analysis basically follows Ref. Q , which re¬ 
views in detail the maximum entropy method applied to 
data of lattice QCD simulation. Here we briefly summa¬ 
rize just several formulae necessary for the later descrip¬ 
tion of our analysis. 

To obtain A{u)) from C{t) by solving the inverse prob¬ 
lem Eq. the maximum entropy method maximizes 
a functional 


p-u}t I 

■ (27) 

This is the continuum type kernel, and an alternative 
form of the kernel which explicitly incorporate the lat¬ 
tice structure was also applied in the literatures [ill^ 
We have seen no signihcant difference in applying the two 
kernels, therefore we shall not further discuss this point. 
In practice, the above integration over the frequency u) is 
replaced by a summation with sufficiently fine discretiza¬ 
tion Aw and a cut off at some maximum value tOmax- 

As already mentioned, one of the main analysis pro¬ 
cedures is the maximum entropy method (MEM) 
Before applying it to finite temperature problems, one 
should verify its applicability under the condition which 
one encounters at T > 0. For this we test how the result 
changes when varying the number of data points used in 
analysis for the correlators at T = 0. The condition to be 
satisfied is that the ground state peak - which we know to 
exist at T = 0 ~ is correctly reproduced, at least at a qual¬ 
itative level. We stress that without such verification, 
the result of MEM may contain uncontrolled artifacts at 
finite temperature. Concerning quantitative questions, 
we find that with reduced number of data points, MEM 
hardly gives a result beyond the qualitative level, for ex¬ 
ample for the width of a peak. 

Another technique in our hand is the standard 
assuming a certain ansatz for the spectral function. This 
method has been applied to a problem of glueballs at fi¬ 
nite temperature . A disadvantage of this approach is 
that one needs to assume a form of the spectral function, 
which introduces a bias. For this purpose, the result of 
MEM can be a good guide. Once a specific ansatz is used, 
the fit gives much more reliable results than MEM for 
the parameters, both in statistical as well as systematic 
sense. 

Therefore MEM and the x^ fit are complementary to 
each other at this stage of computational power, and in 


Q{A]a)=aS[A]-L[A]. (28) 

L[A\ is the likelihood function, 

tiM 

(29) 

where CA{t) is the resulting correlator 12611 for the trial 
A(oj), and V{ti,tj) is the covariance matrix of the data 


1 

= A(iV-l) ^ 

(30) 

with Ck{t) the fc-th sample of the correlator. The stan¬ 
dard x^ fit minimizes this functional L. The Shannon- 
Jaynes entropy 5" [A] is defined as 



A(w) 


to(w) - A(w) log 


(31) 

The function to(w) is called the default model, and should 
be given as a plausible form of A(w). The parameter a 
can be integrated out at the last stage of calculation. 
Following Ref. Q, we use a form 


m(w) = ■ 


(32) 


In the case of point correlators, a natural choice of moM 
is determined according to the asymptotic behavior of 
the meson correlators at large uj in perturbation theory. 
Although such an asymptotic behavior is not observed 
in practical simulation because of the finite lattice cut¬ 
off, this form has been successfully applied to problems 
at zero temperature im For the smeared correlators, 
the situation is more subtle, since the high frequency part 
of the correlator is suppressed by the smearing function. 
We therefore adopt a practical choice: we normalize the 
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smeared correlator so as to produce the same overlap 
with the ground state as the point correlator. Then the 
same normalization is also applied to the correlators at fi¬ 
nite temperatures, and we observe how the result changes 
with the change of muM- 

In the maximization step of Q{A] a) the singular value 
decomposition of the kernel K{t,u!) is used. Then the 
spectral function is represented as a linear combination 
of the eigenfunctions of K{t,uj). The number of degrees 
of freedom of Aiuj) is accordingly reduced to the number 
of data points of the correlator. Although the coefficients 
of the linear combination for A(w) can in principle be 
determined uniquely from the data without introducing 
an entropy term, the small eigenvalues of K(t, oj), 9i, lead 
for the spectral function to a singular behavior; hence 
truncation at some i is practically required M- In MEM, 
the addition of the entropy term stabilizes the problem 
and guarantees an unique solution for the coefficients of 
the eigenfunctions In our analysis, we use only the 
basis functions for which 9i > 10“^^ x 9i holds, where 
9i is i-th largest singular value. This restriction has no 
significant effect on the result. 


C. fit with ansatz for the spectral function 

The standard fit method minimizes the likelihood 
function L, Eq. with an assumed form for the spec¬ 
tral function. A most simple form for fit function is the 
delta function: 


Apoie{Lo;r,m) = rS{u! - m). (33) 



(B 


FIG. 2: The result of MEM analysis for the spectral function 
from the smeared PS correlator composed of free quarks at 
T = 0. 


1-BW form: if the contribution from the high frequency 
part of the spectral function is negligible, a collec¬ 
tive excitation at low energy is expected to be well 
represented by a single Breit-Wigner type ansatz: 

A^oj) = Abw{(^; To, mo, 7o). (36) 

This form is also a good representation for the spec¬ 
tral function of the smeared correlator composed of 
free quarks. 


This form is referred to as pole ansatz in the following. 
At T = 0, a sum of several pole terms should describe 
well the correlators. For the correlators below where 
narrow thermal widths are expected, the multi-pole form 
is also convenient to test this assumption. 

As noted in Sec. m to describe a peak structure with 
finite width we adopt the relativistic Breit-Wigner form 
(referred to as BW form). 


ABw{i^;r,m,j) 


(u;2 - to2)2 -p m272 ’ 


(34) 


This is the same form as Eq. m, with slightly modified 
notation for convenience. This ansatz neglects the uj de¬ 
pendence of r, and hence is valid only for the case that 
the width of the peak, 7, is sufficiently small compared 
with the change of the smearing function in the region of 
interest. 

Combining these two forms, we fit the numerical data 
of correlators to the following ansdtze: 


BW+pole form: although the lowest peak is well repre¬ 
sented with the Breit-Wigner type function, there 
may exist contribution from high frequency region. 
Since the smearing of operator suppress such con¬ 
tribution, remaining effects of this region may be 
represented as a single pole-like term. This is the 
reason that we fit the data to the form 


A{u;) = ABvv(w;ro,TOo,7o) + Apoie{uj-,ri,mi). (37) 

This is the most general ansatz for fit in this 
paper. 

These ansdtze form are the basis for the analysis of the 
lowest peak structure, which corresponds to the ground 
state at T < Tc- The structure in the high frequency 
region is out of the scope of this paper. 


D. Analysis of correlators composed of free quarks 


2-pole form: this form is suitable for description of cor¬ 
relators at T = 0. It is also expected to be a good 
representation of correlators at T < T^, 

A{uj) = Apoieiuj; ro, mo) -k Apoie{uj; ri, mi). (35) 


To learn how to distinguish physical features of the 
correlators from artifacts due to the smearing we apply 
MEM and the fit analyses to the correlators composed 
of free quarks on our lattice, which were discussed in 
sec. imn 
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FIG. 3: The result of fit analysis for the spectral functions 
from the PS correlators composed of free quarks at Nt = 26. 


For the MEM analysis, the fluctuations of the correla¬ 
tors are given by hand in the same manner as the mock 
data analyses in Refs. nm The deviations of the cor¬ 
relators are less than those of the data from the Monte 
Carlo simulation. Figure |21 shows the result of MEM ap¬ 
plied to the smeared correlator in pseudoscalar channel 
composed of free quarks at T = 0. The upper bound 
of t region used for the analysis, t^ax, is varied to ob¬ 
serve how the result of MEM depends on the t-region 
used. With decreasing tmax, the resolution of the spec¬ 
tral function becomes worse, and for tmax £ 16 the spec¬ 
tral function extracted with MEM displays just one peak. 
Therefore for the circumstances specific at T > 0 MEM 
does not have sufficient resolution to distinguish individ¬ 
ual states if they are closer than about 0.05a“^. If the 
correlator is composed of almost free quarks, the width 
extracted with MEM may be of order of 0.05-0.1 in tem¬ 
poral lattice units. 

We also apply the m analysis to the correlators 
composed of free quarks. In this case the errors of the 


same size as in Monte Carlo simulation are just put on 
the correlators without fluctuating them. Figure 13 shows 
the result of the fit of the correlators composed of 
free quarks at Nt = 26 using the 1-BW ansatz. The 
imin dependence of mass and width (with fixed tmax = 
13) indicates that the single BW ansatz describes rather 
well the correlators. As is observed in Figure El if the 
correlator is composed of free quarks, the x^ fit gives 
sizable difference in mass and width parameters for the 
smeared and half-smeared correlators. This dependence 
is in agreement with the fact that the propagator should 
only show the two free quark cut and no particle-like 
excitations and indicates that by testing the dependence 
of the result on the smearing function we can distinguish 
physical effects from artifacts due to smearing. 


V. SETUP OF NUMERICAL SIMULATION 
A. Lattice setup 

The zero temperature lattice used in this paper is 
the third one of Ref. 1^: a quenched lattice of size 
20^ X 160, generated with the standard plaquette action 
with (/3,7g) = (6.10,3.2108). These coupling and bare 
anisotropy correspond to the renormalized anisotropy 
^ = aajar = 4 within 1% accuracy [^, and the spatial 
lattice cutoff a~^ = 2.030(13) GeV set by the hadronic 
radius ro [l^. At T = 0, 500 configurations are gener¬ 
ated with the pseudo-heat-bath update afgorithm, each 
separated by 2000 sweeps after 20000 sweeps for ther- 
malization. The mean-field values are defined as the av¬ 
erage values of link variables in the Landau gauge, and 
obtained as Ua = 0.8059(1) and Ur = 0.9901. 

To determine the critical temperature, we measure the 
Polyakov loop susceptibility at Nt = 27, 28, and 29 at 
f3 = 6.10, and in addition, at several values of [3 (with 
corresponding values of xg) around (3 = 6.10 at fixed 
Nt = 28. At /3 = 6.05 the lattice scale set by ro is a~^ = 
1.892(10) GeV, which together with a~^ at (3 = 6.10 
determine the scales at the other values of [3 by linear 
interpolation. The susceptibility peaks at about [3 = 
6.10 and Nt = 28. The critical temperature is obtained 
as Tc ~ 290 MeV with 10 MeV of roughly estimated 
uncertainty. 

The charmonium correlators at T > 0 are measured 
for two values of temporal lattice extent, Nt = 32 and 
26. Corresponding temperatures are 0.88rc for Nt = 32, 
and 1.08rc for Nt = 26. For brevity, these temperatures 
are hereafter referred to as 0.9Tc and l.lTc, respectively. 
Thus the temperatures treated in this paper are in the 
vicinity of the transition. At each of these two Nt’s, 
we generate 1000 configurations each separated by 500 
pseudo-heat-bath sweeps after 20000 sweeps for thermal- 
ization. 

At Nt = 26, we find that the configurations almost 
stay in a single Polyakov loop sector during the whole 
updating history. Transitions to other sectors occurred 
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TABLE I: The spectrum at zero temperature determined 
from the point and the smeared correlators. The masses in 
PS and V channels are determined with the two-pole ht. 


Correlator 

State 

mps 

mv 

Point 

ground 

0.36835(37) 

0.37748(49) 


first exc. 

0.449(31) 

0.463(42) 


fit range 

45-80 

50-80 

Smeared 

ground 

0.36856(9) 

0.37769(12) 


first exc. 

0.500(22) 

0.479(23) 


fit range 

17-80 

15-80 


only once after more than 470k sweeps. In Ref. Q it 
was reported that the light mesonic correlators behave 
differently on configurations in different sectors. We re¬ 
gard quenched lattices as approximations to lattices with 
dynamical quarks, on which Z 3 center symmetry is ex¬ 
plicitly broken and Polyakov loop prefers to stay on the 
real axis. Therefore we transform all the configurations 
to the the real sector of the Polyakov loop. 


B. Charmonium correlators 

As mentioned in Sec. ITTT Rl we use (k, 7 ^) = 

(0.112,4.00), which correspond to mg ~ 0.98 GeV, 
roughly the charm quark mass. The statistical errors 
for the results are estimated with the jackknife method 
with appropriate binning. 

In the next section, we start our analysis at T = 0 with 
the examination of the point correlators. On each con¬ 
figuration, we calculate the point correlators four times 
with four different source points: t = 0 and 80 at spatial 
site X = (0, 0,0) and t = 40 and 160 with x = (10,10,10). 
Then these correlators are averaged with appropriate 
shifts in t-direction. At T = 0 the standard fit of data 
to a multipole form works well. To fix the general picture 
we anticipate on the analysis of the next sections and list 
in Tabled the result of two-pole fits of the point correla¬ 
tors in PS and V channels. The fit ranges [tmin, tmax] are 
determined by varying tmin with fixed tmax = 80 and ob¬ 
serving the stability of the fit parameters. The results of 
the ground state masses are very close to the experimen¬ 
tally observed charmonium masses. On the other hand, 
the hyperfine splitting is with about 74 MeV smaller than 
the experimental value, mj/^ — niri^ ~ 117 MeV. This is 
a well-known feature of the quenched lattice simulations 
and is considered mainly due to quenching, although lat¬ 
tice artifacts in the charmonium system can also play an 
important role |3f| . 

The smeared correlators are the main target of this 
paper for reasons we already explained. The smearing 
parameters are described in Sec. mm The result of a fit 
to 2-pole form is also shown in Tabled It is well-known 
that the correlators with smeared sink are quite noisy. 
To reduce the noise, we measure sixteen correlators with 



t 

FIG. 4: The effective masses for the smeared correlators in 
PS and V channels at three temperatures. 


different source points on each configuration, and average 
them. We change the spatial center of the smearing func¬ 
tion, as well as the time slice, to reduce the correlation 
as much as possible. In contrast to a naive expectation, 
we find this procedure efficient. In fact, at T = 0, the 
statistical errors are reduced by a factor of about 0.3 in 
the region t = 8-16, which is the most important region 
for the present analyses, in both the PS and V channels. 
This reduction of errors corresponds to an increase in 
statistics of about 10 times. Almost the same amount of 
reduction of errors is observed at Nt = 32. At W = 26, 
the errors are reduced by factors of 0.40-0.45 in the range 
t = 8-13. This way of reducing statistical errors will be 
advantageous in the case of dynamical simulations. 

For the present kind of analysis, whether the analysis 
is efficient or not may significantly depends on the statis¬ 
tics. In the temporal region t = 8-16, the statistical er¬ 
rors of point correlators at T = 0 are about 0.12-0.18% 
and 0.08-0.15% for PS and V channels, respectively. In 
the same region, the smeared correlators have statistical 
errors of 0.14-0.17% and 0.15-0.20%. Therefore at T = 0 
the errors are of the same order for the point and smeared 
correlators. At finite temperature, with 1000 configura¬ 
tions, the statistical errors are of roughly the same size 
as at T = 0. We also note that the correlation between 
the correlators at neighboring time slices is stronger for 
the point correlators than for the smeared ones. 

To obtain an impression of the physics to be expected 
at various temperatures we show in Figure 01 the effec¬ 
tive mass plot for the smeared correlators, with meffit) 
defined through 

Cjt) ^ cosh [meffit)iNt/2-t)] 

Cit^-l) cosh [TOe//(t)(Vt/2 - t - 1)]' 

If the correlator is dominated by a single (stable) state, 
meffit) shows a plateau. In Fig. ^ the effective masses 
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at T = 0 show plateaus beyond t ~ 16. The effective 
masses at Nt = 32 (T ~ 0.9Tc) show almost the same 
behaviors as at T = 0. This implies that the charmo- 
nium states are almost unchanged at this temperature, 
compared to those at T = 0. In contrast, me//(t)’s at 
Nt = 26 (T ~ l-lTc) show quite different behavior: they 
decrease rapidly compared to those at T < Tc, and in the 
large t region {t > 10 ) the pseudoscalar and vector chan¬ 
nels are almost degenerate. Therefore a qualitatively dif¬ 
ferent behavior of the spectral function is expected below 
and above T^. These features are consistent with obser¬ 
vations in earlier work Q. 

VI. ANALYSIS AT ZERO TEMPERATURE 
A. Result of MEM for point correlators 

We start with MEM analysis of the point correlators 
at zero temperature. 

First we discuss the default model function. Follow¬ 
ing previous applications of MEM to lattice data nm, 
we adopt as moM in Eq. 1221 ) the asymptotic values of 
correlators in the perturbative QCD, 

» Aqcd) = ^(l + r2^), (39) 

where ri =3/2 and r 2 = 11/3 for PS channel, and ri = 1 
and r 2 = 1 for V channel. As value of the strong coupling 
constant Us, we use Us = 0.2 at 8 GeV from Ref. as 
a typical value at the temporal cutoff of the present lat¬ 
tice. The matching of the lattice theory with the contin¬ 
uum theory is performed through the tadpole-improved 
tree-level, and hence the renormalization of the quark 
field is represented only by the KLM normalization fac¬ 
tor (EH). These settings give for the parameter of the 
default model function uidm = 4.1 for PS channel, and 
2.4 for V channel. If MEM works as a method to extract 
the spectral function without assuming a specific form, 
the result should not be sensitive to the choice of this pa¬ 
rameter. In addition, the observed asymptotic behavior 
of the spectral function from lattice data around the cut¬ 
off is actually different from that of Eq. El- We there¬ 
fore observe the dependence of the result on the value 
of moM by changing ttidm by factors of 10 and 0 . 1 , to 
verify the insensitivity of the result to muM- 

The further setting of MEM is as follows. In all 
cases, the correlator at t = 0 is not used for the anal¬ 
ysis. The minimum of t used depends on the type of 
analysis, in most cases tmm = 1 is adopted. The fre¬ 
quency u is discretized with Aw = 0 . 002 , in the region 
= [0.001,4.0]. The dependence of the result 
on the parameters for w is sufficiently small around the 
adopted values. 

We first examine how the restriction of degrees of free¬ 
dom of the correlator affects the reconstructed spectral 
function. Three parameters tmm, tmax and tsep are in¬ 
troduced, that is MEM is applied for the correlators at 
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FIG. 5: The result of MEM analysis of the point correlator in 
PS channel at T = 0. The figure shows the tmax dependence 
of the resultant spectral function. 



FIG. 6: The result of MEM analysis of the point correlator in 
PS channel at T = 0. This figure shows the tsep dependence 
of the result of the spectral function with fixed tmax = 48. 

times tmin < t < tmax in every tsep time slices. We ex¬ 
amine two types of restrictions of the numbers of degrees 
of freedom: 

(!)■ tmax is varied as 48, 24, 16, and 12, with fixed tmin — 
1 and tsep = 1. The corresponding numbers of data 
points are the same as tmax- 

(II). tsep is varied as 1, 2, 4, and 8, with fixed tmax = 48. 
tmin is varied accordingly as tmin = tsep- Corre¬ 
spondingly, the numbers of data points are 48, 24, 
12 , and 6, respectively. 

From the point of view of the basis functions of the sin¬ 
gular value decomposition in the spectral function space. 
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case (I) reduces the number of basis functions while keep¬ 
ing the functions unchanged. On the other hand, case 
(II) dilate the basis functions by a factor of tsep simul¬ 
taneously reducing their number. Reliability of the re¬ 
sult at finite temperature requires the obtained spectral 
function to be stable under the type (I) restriction. Al¬ 
though larger values for tmax are possible, tmax = 48 
is sufficiently large and has appropriate physical range 
{tmaxCiT — 1-2 fm) for the present purpose. 

The result is shown in Figs. 0 and ini for the pseu¬ 
doscalar channel. For the vector correlator a similar re¬ 
sult is obtained. In the case of tmax = 48 and tsep = !> 
i.e. without restriction, the reconstruction of the spectral 
functions seems successful, and their fundamental fea¬ 
tures appear to be the same as in previous works nm 
The peak positions of the ground and first excited states 
are consistent with the result of fit summarized in Ta¬ 
ble ^ Here we do not discuss the higher excited states 
because of uncertainties not only in MEM but also in the 
fit analysis. 

Figure El displays the result under restriction condi¬ 
tions (I). Decreasing tmax, the spectral function becomes 
increasingly different from that with tmax = 48. For 
tmax < 24, the first excited state peak does not appear 
at the correct place, and for tmax 16 even the ground 
state peak is located at an incorrect place and has broad 
width. These lowest peaks for different tmax’s are all sig¬ 
nificant in the sense of the error analysis of MEM Q. 
Therefore tmax < 16 is not acceptable for extracting reli¬ 
able information of the ground state from the point cor¬ 
relators. Since at finite temperature we are restricted to 
tmax < 16 because of short temporal extent this is the 
reason why we abandon using the point correlators and 
apply the smearing technique. 

Eigure E) shows the result with restriction condition 
(II). In this case, the reconstructed spectral functions 
are rather stable, while the sharpness of peaks is lost 
by increasing tgep- Even with 6 degrees of freedom (as 
a linear combination of 6 functions), the reconstructed 
spectral function at least exhibits the ground state peak 
at the correct position. To summarize, if one has a tem¬ 
poral region of the correlator of the order of 1 fm, the 
spectral function can be rather nicely reconstructed even 
with less than 10 degrees of freedom. 

From these observation, we conclude that with 0(10) 
degrees of freedom the crucial condition is determined 
by the region of t where the correlators are measured. 
tmaxCir of the Order of 1 fm seems necessary for a reliable 
extraction of the spectral function from the point corre¬ 
lator. Increasing the number of degrees of freedom above 
10 with fixed physical range of [tmim, tmax] improves the 
situation only slightly. 

Now we consider the effect of the statistics. We apply 
the same analysis of the case (I) restriction to the correla¬ 
tors averaged over first 100 configurations. The observed 
features of tmax dependence are essentially the same as 
with 500 configurations. Therefore at the present level of 
statistics (order of several hundreds), increasing statis¬ 


tics does not seem to improve drastically the situation 
of the tmax dependence. We also compare the correla¬ 
tors averaged over 100 configurations randomly selected 
from the total of 500 configurations. A comparison of 
the results with tmax = 48 from five such correlators 
shows that the first excited state does not always ap¬ 
pear at the same place as with 500 configurations, but 
deviates in a range of w = 0.44-0.52. Thus with less 
statistics only the ground state peak can be trusted even 
with tmax = 48. These observations indicate the required 
statistics for each situation. 

Finally we briefly discuss the dependence on the de¬ 
fault model. We compare the results obtained with ttidm 
in Eq. (E3l multiplied by 10 and 0.1 with otherwise fixed 
parameters, tmax = 48 and tsep = 1- These results of the 
spectral function display the peaks at the same places as 
with original itidm, while increasing moM the peaks be¬ 
come sharper. Although present analysis is quite simple, 
we conclude that the default model dependence of the re¬ 
sult of MEM is not severe, at least for a rough estimate 
of the shape of the spectral function. 

The most important conclusion in the analysis of the 
point correlator is that we need the order of 1 fm for the 
range of t where the correlators are measured, for reliable 
extraction of the spectral function from the correlator 
with present level of statistics. This requirement cannot 
be fulfilled at finite temperature, and a brute force appli¬ 
cation of MEM to the point correlators at T > 0 would 
produce an unreliable result. This result warns us against 
physical significance of the results in Refs. p^ftllTlIlMIl^ . 
This difficulty may originate in that the point correla¬ 
tor at short distance contains the contribution from wide 
range of frequency up to the lattice cutoff. We there¬ 
fore give up to analyze further the point correlators, and 
concentrate our attention on the smeared correlators. 


B. Result of MEM for smeared correlators 

We now turn to MEM analysis of the smeared corre¬ 
lators. Since the smearing technique reduces the high 
frequency part, there is no a priori choice for the form of 
the default model. As described in Sec. nvH we use the 
same default model function as for the point correlator, 
with the normalization of the smeared correlator which 
gives the same overlap with the ground state as the point 
correlator. The default model dependence is examined in 
the same manner as for the point correlators. Other pa¬ 
rameters are the same as for the point correlator, except 
that ujmax = 2.0 is adopted, since the high frequency 
region is sufficiently suppressed. 

Eirst we examine the dependence of the result on the 
numbers of degrees of freedom, using the same restric¬ 
tion conditions as for the point correlators. Since the 
smearing technique enhances the low energy part, we only 
consider whether the ground state peak is correctly re¬ 
produced or not. The result for the restriction condition 
(I) is displayed in Fig. [Tj In the both cases of PS and 
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FIG. 7: The spectral functions at Nt = 160 determined with 
MEM. The top and bottom panels show the tmax dependences 
of obtained spectral functions in PS and V channels, respec¬ 
tively. 


V channels the place of the ground state peak is stable 
with decreasing tmax, although the width is gradually 
broadened. Since the ground state peak evidently has no 
physical width at T = 0, this hctitious width is a system¬ 
atic uncertainty of MEM due to the insufficient number 
of basis functions in the spectral function space. There¬ 
fore, in particular at finite temperature, an estimate of 
the physical width requires a careful analysis combined 
with other procedures. The instability of the first excited 
state peak under changing tmax is explained by the small¬ 
ness of the contribution to the correlators from energies 
above the ground state. In fact, the two-pole fit analysis 
indicates that the overlap of the first excited state with 
the smeared correlator is about 7% (6%) of that of the 
ground state for the PS (V) channel. Thus the instabil¬ 
ity of the excited state peaks has no signihcance for the 
following analysis. We conclude that for the smeared cor¬ 



FIG. 8: The default model dependence of the extracted spec¬ 
tral function with MEM from the smeared correlators in PS 
channel. The figure shows the result of MEM with three val¬ 
ues of irtDM, the parameter of the default model function 
Eq. 18211 . 


relators, MEM with tmax — 12 would work for a rough 
estimate of the structure of the spectral function in the 
low energy part, such as examination of whether strong 
ground state peak exist, while a quantitative evaluation 
of width, for example, is difficult. In the case of the 
restriction condition (II), a similar tendency as for the 
point correlators is observed: the position of the ground 
state peak is stable with increasing tgep, while the sharp¬ 
ness of the peak is gradually lost. 

The rriDM dependence of the result is examined in two 
cases of tmax, tmax = 48 and 12, by using rescaling factors 
10 and 0.1. In the case of tmax = 48, the ground state 
peak is strong, and changing uidm does not cause a large 
effect other than the decrease of the peak height with 
slight increase of the width when decreasing to_dm- Fig-IHl 
shows the case of tmax = 12 for the PS channel. Although 
the peak shape is broadened with decreasing muM, the 
essential features of the result are stable. A similar result 
is observed for the V channel. 

Let us summarize the MEM analysis of the smeared 
correlators at T = 0. In contrast to the case of the point 
correlators, MEM with restricted numbers of degrees of 
freedom works up to tmax = 12, i.e. tmax^r — 0.3 fm, 
at least for a rough estimate of the shape of the spec¬ 
tral function in the low energy part. For this purpose, 
systematic uncertainties in the default model parameter 
do not appear significant. This is an encouraging result 
for an application of MEM to finite temperature as a 
precedent analysis to the Ft- 
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FIG. 9: The result of the x fif analysis at Nt = 160. The 
top and bottom panels show the results for mass and width 
of the ground state peak, respectively. 
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FIG. 10: Spectral functions at Nt = 32 determined with 
MEM. 


C. Result of fit analysis 

The results with 2-pole ansatz have been given in Ta¬ 
ble Q] For completeness, we also apply the fits to the 
smeared correlators with forms other than the 2-pole fit. 
The result is shown in Fig. 1^1 As fit range [tmimtmax] 
we fix tmax = 13, considering the severest case at T > 0, 
and observe the tmin dependence of the result. The top 
and bottom panels of Fig.|^show the results of mass and 
width of the ground state peak, respectively. 

As is expected, the 2-pole form describes well the cor¬ 
relators, and the other two fit forms are consistent with 
the 2-pole fit. In the cases of 1-BW and BW-Ppole hts 
the tmin dependences of the widths are signihcant. In 
the bottom panel the results for the width decreases as 
tmin increases, and seem to approach zero. Therefore no 
indication of a finite width is found at T = 0, as it should 
be. 


VII. ANALYSIS AT r < Tc 

Now the two analysis procedures are applied to the 
correlators at T ~ 0.9Tc {Nt = 32). In this and the next 
sections we no longer discuss the point correlators and 
focus only on the smeared ones. 


A. Result of MEM analysis 

We begin with the MEM analysis, in accordance with 
our strategy. The parameters concerning MEM are al¬ 
most the same as for the smeared correlators at T = 0. 
The range of t is fixed to tmin = Ij tgep = 1, and 
t — 

^max — -L'-'* 

The result is displayed in Fig. cni for three values of 
'm-DM- For the reasons already mentioned we focus on 
the ground state peak and do not discuss in detail the 
high frequency part of the spectral function. As appar- 
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ent in the figure a ground state peak appears in both the 
PS and V channels. Although with decreasing toom the 
widths of the ground state peaks increase, the positions of 
the peaks are stable and almost the same as at T = 0. In 
the present case, since there is no intrinsically advisable 
value for moM beyond an order estimate, this ambigu¬ 
ity of the width of the peak should be considered as an 
uncertainty of MEM applied to the smeared correlators. 
We also perform the same analysis with less statistics, 
500 configurations, to see how this result depends on the 
statistics. The result is essentially the same, and thus 
statistically stable. 

These results support the assumption that the mesonic 
ground states are persistent up to this temperature, with 
almost the same masses as at T = 0. The width of 
the ground state peaks is small while finite, although 
it strongly depends on the value of rriDM- Whether the 
width is physically finite or not should therefore be ex¬ 
amined with the fit analysis. 


B. Result of fit analysis 


According to the result of MEM we can assume that 
the spectral function of the smeared correlator at this 
temperature has a peak structure at almost the same 
mass as at T = 0 and with small width. Therefore it is 
expected that the fits with the three types of fit forms, 
2-pole, 1-BW, and BW-|-pole forms can clarify the low 
energy structure of the spectral function. 

The results are summarized in Fig. im The top panel 
displays the dependence of the mass parameters for the 
ground state peak on the lower bound of the fit range, 
tmin- The upper bound is fixed to t^ax = 16. It is appar¬ 
ent that the fit to the 2-pole form exhibits stable results 
both for the PS and V channels beyond tmin = 3. The 
result of fit to 1-BW form exhibits no plateau, indicating 
that this form does not explain the whole structure of the 
correlators. However, the values seem to approach to the 
corresponding results of the 2-pole fits. The result of fit 
to the BW-|-pole form is consistent with those of 2-pole 
fit at tmin ^ 2. 

The bottom panel of Fig. ^2 shows the result for the 
width parameter. In the case of the fit to the 1-BW form 
the value of width gradually decreases with increasing 
tmin ■ This behavior is consistent with a vanishing width. 
The fit to BW-|-pole form also indicate that the width is 
consistent with zero. Therefore, there is no indication of a 
finite width for the spectral function at this temperature. 

All the results of fits to the three forms indicate that 
the ground state peak is well described by a strong peak 
with vanishing width, i.e. a pole-like structure, and the 
associated mass is almost the same as at zero tempera¬ 
ture. 
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FIG. 11: The result of x^ fit analysis at Nt = 32. The top 
and bottom panels show the results of the mass and width 
parameters for the ground state peak, respectively. 


VIII. ANALYSIS AT T > 

A. Result of MEM analysis 

In this section, we analyze the correlators at T ~ l.lTc 
{Nt = 26). 

The setup of MEM analysis is the same as at T ~ 
0.97)., except for tmax = 13. Figure d shows the re¬ 
sult. In both the PS and V channels there appear strong 
peak structures around lo ~ 0.4. This peak, hereafter 
called “the ground state peak” for simplicity, appears at 
almost the same position as at T ~ 0.9 but with larger 
width. Although, as a common tendency, the peak be¬ 
comes sharper as the default model parameter ttidm in¬ 
creases, the position of the peak is unchanged. It is also 
verified that the result is essentially the same with less 
statistics, namely 500 configurations. Therefore it is pre¬ 
sumable that the spectral function still has a peak struc¬ 
ture at almost the same position as below Pc- 
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FIG. 12: The spectral function at Nt = 26 determined with 
MEM (smeared correlators). 


This result should be compared with the case where 
the correlators are composed of free quarks, considered 
in Sec. imn In the latter case, as shown in Sec. 

MEM with the present number of degrees of freedom may 
produce the same feature of the spectral function but 
with a lower peak position. The above result therefore 
does not exclude the possibility that the correlators are 
composed of two weakly interacting quarks with rather 
large effective mass. 

To judge between almost free quarks and genuine 
bound-state-like structure we analyze the half-smeared 
correlator. We repeat the same analysis as applied to 
the smeared correlators for the half-smeared correlators 
and verify that MEM works up to tmax = 12 on almost 
the same level as for the smeared ones. As discussed 
in Sec. IIVT^ if the correlator is composed of two almost 
free quarks the extracted spectral function should exhibit 
now a wider peak at a higher position. Figure El shows 



00 


FIG. 13: The results of MEM for the spectral function of the 
half-smeared correlators at Nt = 26. 


the result of MEM analysis for the half-smeared correla¬ 
tors. The analysis is performed in the same manner as 
for the smeared correlators. Comparing with the result 
for the smeared correlator, the peak position for the half- 
smeared one is almost unshifted, while the width of the 
peak tends to broaden slightly. The latter effect , how¬ 
ever, can also be easily explained as an effect of default 
model uncertainty. This result contradicts therefore the 
assumption that the peak structure is an artifact of the 
smearing. Hence the MEM analysis supports a bound- 
state-like structure at this temperature. 


B. Result of fit analysis 

According to the result of MEM, it is reasonable to 
assume that the spectral function at this temperature 
also have a peak structure, similarly to the T < 
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FIG. 14: The result of fit analysis for the smeared corre- FIG. 15: The result of x fit analysis for the half-smeared 


lators at Nt = 26. 


correlators at Nt = 26. 


case. Therefore the three forms for the x^ fit analysis 
are still reasonable assumptions to investigate the low 
energy structures of the correlators. 

Fjgure im shows the result of the x^ fit analysis for the 
smeared correlators. The top panel shows the tmin de¬ 
pendence of the mass parameters. In contrast to the case 
below Tc, the masses from 2-pole fit no longer approach 
those of the other two fit forms but fall monotonously. 
The masses from 1-BW and BW-|-pole fits exhibit consis¬ 
tent behaviors, indicating that these fits represent better 
the correlators. Actually the values of x^/^df are con¬ 
sistently fluctuating around unity at tmin > 8 for 1-BW 
fit, and in the whole range of tmin for BW-|-pole fit. In the 
BW-|-pole fit, the mass parameters for the second peaks 
take values around 0.9 and 1.0 for PS and V channels, re¬ 
spectively, which is roughly consistent with the result of 
MEM. The consistency of 1-BW and BW-Lpole fits also 
holds for the width parameters, as displayed in the bot¬ 
tom panel of Figure ITU These results indicate therefore 
that the widths associated with the ground state peaks 


are finite for both the correlators in PS and V channels. 

We now analyze the half-smeared correlators to exam¬ 
ine whether this observation of peak structure is an ar¬ 
tifact of the smearing. The analysis is performed in the 
same manner as for the (fully) smeared correlators. The 
result is displayed in Figure IT^ The mass parameters for 
the half-smeared correlators (top panel) show a similar 
tendency as for the smeared correlators, although the re¬ 
sult of 1-BW and BW-Lpole fits approach slightly larger 
values than those of smeared correlators. The shifts of 
the masses are at most 5% and can be explained with the 
frequency dependence of the overlap R{lu) in Eq. (H^ . A 
similar tendency is also found for the width parameters in 
the bottom panel of Fig. lTKl which appear even more con¬ 
sistent with those of smeared correlators than the mass 
parameters. Therefore we find no significant shift of the 
peak position or broadening of the peak width, such as 
expected for the case with almost free quarks discussed 
in Sec. lIIK)! 

Figure ITTI illustratively compares these spectral func- 
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FIG. 16: The spectral functions obtained with fit analysis 
at Nt =26. For the correlators composed of free quarks only 
the result of BW fit is displayed, while a BW+pole fit is used 
for the correlators measured in the simulation. 


tions extracted with the fits from four kinds of corre¬ 
lators: the smeared and half-smeared correlators, numer¬ 
ically measured and composed of free quarks. The results 
of the numerical simulation are of the BW-l-pole fit with 
tmin = 4. For the free quark case the results of 1-BW 
fit with tmin = 8 is displayed. We note that the normal¬ 
izations of the spectral functions are not significant in 
the present analysis. In both the PS and V channels the 
spectral functions for the numerically observed correla¬ 
tors strongly peak near lo = 0.4 with small widths. The 
effect of the smearing seems not significant. In contrast, 
the spectral functions for the correlators composed of free 
quarks are strongly dependent on the applied smearing. 
Both a sizable shift of mass and a broadening of width 
are observed in this case. It is difficult to exclude com¬ 
pletely the possibility that the correlators are composed 
of almost free quarks with rather large effective mass 
and nontrivial dispersion relation. However, it is quite 


suggestive that a nontrivial structure indicating the ex¬ 
istence of bound-state-like structure in the low energy 
region of spectral function subsists at this value of T. 

Summarizing the results of MEM and fit analysis we 
conclude that the charmonium correlators in PS and V 
channels at T ~ l.lTc possess a nontrivial peak structure 
in the low energy region. As representative values of toq 
and 7 o, we quote the results of the BW-bpole fit analysis 
with tmin = 4 for the smeared correlators: 

PS : mo = 3.06(2) GeV, 70 = 0.12(3) GeV, (40) 
V : mo = 3.13(2) GeV, 70 = 0.21(2) GeV. (41) 

The quoted errors in mo and 70 are only the statistical 
ones and do not include the error in a~^. As obvious 
from the discussions in this sections, these results should 
contain systematic uncertainties of the order of 5% due 
to the analysis procedures, apart from other uncertain¬ 
ties such as finite lattice artifacts and quenching effects. 
Compared with the result at T ~ 0.9Tc, the widths at 
T ~ l.lTc is sizable, indicating a genuine temperature 
effect in the deconfined phase. On the other hand, the 
masses are almost unchanged. This result is contradic¬ 
tory to the absence of bound states argued from the po¬ 
tential model approach | 2 , • 


IX. CONCLUSION 

The main goals of this paper were (1) to elucidate the 
technical problems in extraction of the spectral function 
from lattice data of t-correlators, and ( 2 ) to investigate 
the temperature effect on the spectral function near the 
deconfining phase transition. In the following, we sum¬ 
marize and discuss the results obtained in this paper. 

As techniques to extract the information on the spec¬ 
tral function we examined maximum entropy method 
(MEM) and the standard fit assuming suitable forms 
of the spectral function. It is essential to check the relia¬ 
bility of the applied methods to the systems in question. 
Our condition for reliability at finite temperature is that 
the methods reproduce the correct form of the spectral 
function at T = 0 when the t-interval is restricted to the 
one forced on us at T > 0. We examined MEM by ap¬ 
plying it to the point and smeared correlators. We find 
that MEM does not meet this requirement for the former, 
while it does for the latter. This is understandable, since 
the smearing enhances the low energy part of the spec¬ 
tral function, which is what we are interested in, while 
a much wide energy region contributes to the point cor¬ 
relators. We note that whether MEM correctly works or 
not depends primarily on the extension of the physical t- 
region. In particular, for the point correlators a region of 
t of O(lfm) is necessary. For the smeared correlator, this 
condition is much relaxed. Therefore only the smeared 
correlators were analyzed at finite temperatures. 

Since MEM is ambiguous concerning the quantitative 
detail of the spectral function, the latter should be also 
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analyzed with other procedures. As such a procedure, 
the fit is a reasonable candidate, since MEM already 
gives a hint for a suitable ansatz for the spectral function. 
With several assumed forms and examining the fit range 
dependence of the resultant values for the parameters, 
this procedure gives us more quantitative information on 
the properties of spectral function. We emphasize that 
MEM and x^ fit analysis as used here are complementary 
to each other. This two-step approach actually worked 
for analyses of the correlators at finite temperature, as 
well as at T = 0. 

Now we discuss the physical implications of our results 
below and above the critical temperature. We remind 
that our numerical simulation was performed without dy¬ 
namical quark effects. 

Below the deconfining temperature, at T ~ 0.9Tc, the 
reconstructed spectral function has a strong peak corre¬ 
sponding to the ground state, with almost the same mass 
as at T = 0 and narrow width consistent with zero. In 
contrast to the potential model analysis Q, the charmo- 
nium mass is not changed up to this temperature. Similar 
tendencies have been reported in previous lattice QCD 
calculations for the mesonic channels 0, 0, while siz¬ 
able reduction of mass has been observed for glueballs 
Considering the rather quantitative success of the 
potential model approach for the charmonium systems at 
T = 0 [^, it is important to explain this discrepancy. 

At T cz l.lTc we observed an indication that the spec¬ 
tral functions still has strong peaks at almost the same 
positions as T < T^, and with widths of about 0.12 and 
0.21 GeV for PS and V channels, respectively. This re¬ 
sult presumably indicates the existence of quasi-stable 
bound-state-like structures persistent up to this temper¬ 
ature. The possibility of observing correlators composed 


of almost free quarks (but of large effective masses) is, 
however, not completely excluded. Finally, also above 
the observations are not in accord with the expectation 
from the potential model approach |^. This result im¬ 
plies that the plasma phase has a nontrivial structure at 
least near the critical temperature. 

For a more definite understanding of hadronic corre¬ 
lators at T > 0 more studies containing dynamical sim¬ 
ulations are necessary, as well as calculations in wider 
range of temperatures. The techniques adopted in this 
paper should be applicable to these situations. It is also 
important to repeat the same sort of analysis in the light 
hadron sectors, for which a nontrivial structure of corre¬ 
lators above Tc was also reported 0. 
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